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مدل‌های 00۷ به‌طور وسیع برای ارزیابی تغییر اقلیم در یک مقیاس جهانی استفاده می‌شود؛ اما خروجی 
این مدل‌ها برای ارزیابی تغییرات اقلیمی در سطح محلی و منطقه‌ای کافی و دقیق نیست. در این مقاله با 
استفاده از مدل 515۷ خروجی مدل تغییر اقلیم 2 را در منطقه مورد مطالعه به‌وسیله داده‌های 
مشاهداتی ایستگاه تبریز که دارای آمار بلندمدت اقلیمی است. ریزمقیاس نموده و با در نظر گرفتن 
سناریوهای تغییر اقلیم 18۳2.6 30۳4.5 و 80۳8.5 برای دوره‌های آیندهٌ ۲۰۱۰-۲۰۳۹ ۲۰۰-۲۰۹ 
و ۲۰۷۰-۲۰۹۹ تغییر اقلیم منطقه مورد نظر مورد ارزیابی قرار گرفته است. مشاهدات روزانه حداقل و 
حداکثر دما؛ بارش برای دوره پایةٌ ۱۹۱۰-۱۹۹۰ به عنوان ورودی وارد مدل شده است. نتایج حروجی مدل 
ریزمقیاس نشان می‌دهد که در دوره‌های آینده دما در ایستگاه تبریز بر اساس سه سناریوی مورد بررسی 
افوایکن واه بافته این افدایی برع شوه ۲۲ ۲۰۷۰۵۱۵۹ یو تشر وراه بوخ ون 
ایستگاه تبریز به‌طورکلی بارش در سه سناریوی مورد بررسی برای دو دوره ۲۰۱۰-۲۰۳۹ و ۲۰۷۰-۲۰۹۹ 
کاهش و برای دوره ۲۰۶۰-۲۰۹ افزايش می‌یابد. همچنین بارش به‌طورکلی در فصل زمستان افزایش و 
بقیه فصول با کاهش بارش مواجه خواهد بود. تغییرات میانگین حداقل دمای ایستگاه تبریز در کلیه ماه‌ها 
به‌غیراز ماه توامبر و دسامبر دی دوره‌فای, ایثده افزایش داشته است. حداقل دما در سه سناریوی مورد 
بررسی برای سه دوره مورد مطالعه افزايش می‌یابد. همچنین حداقل دما به‌طورکلی در تمام فصول افزایش 
می‌یابد که در فصل تابستان تا ۸ درجه نیز افزایش دما مشاهده می‌گردد. 


کلیدواژه‌ها: بارش دماء تغییر اقلیم نبریز. 


۲۱۵۵11: 2٩21017001 0 ۰٩۹۱۲۲۶۱۹۹۵۸ نویسنده مسئول:‎ ۱ 
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۱. مقدمه 


تغییر اقلیم یکی از گسترش يافته ترین موضوعات مورد بحث در رشته اقلیم‌شناسی در دو دهه اخیر است. 
انسان به عنوان جزئی از سیستم اقلیمی نقش بسزایی در رفتار اقلیم دارد. بخصوص در عصر حاضر با افزایش 
جمعیت و نیاز انسان به منابع آب و غذا و درنتیجه افزایش اراضی کشاورزی از بین رفتن جنگل‌ها 
بیابان‌زایی و افزايش استفاده از انرژی‌های فسیلی باعث دگرگونی در سیستم اقلیمی شده است. طبق چهارمین 
گزارش ارزیابی تغییر اقلیم هیئت بین دول تغییر اقلیم (1۳0/67) پیش‌بینی‌های مدل‌های اقلیم جهانی برای قرن 
بیست و یکم نشان می‌دهند که گرمایش جهانی به شتاب خود ادامه خواهد داد حتی اگر بشر بتواند جلوی 
انتشار گازهای گلخانه‌ای را بگیرد. پیش‌بینی می‌شود که تا سال ۲۱۰۰ تغییرات متوسط دمای جهانی از ۱/۸ 
درجه به ۶ درجه سانتی گراد افزايش یابد و سطح متوسط دریاها بین ۰/۱۸ تا ۰/۵۹ متر بالا بیاید. فراوانی و 
حجم رخدادهای حدی نیز از قبیل سیل. خشکسالی و موج‌های گرمایی با افزايش دمای متوسط جهانی 
گسترش می‌یابد. اثرات تغییر اقلیم بر کل سیاره یکسان نیست. بعضی مناطق زمین حساسیت بیشتری نسبت 
به تغییر اقلیم در مقایسه با مناطق دیگر دارند. به‌عنوان‌مثال مناطقی که دارای اقلیم مدیترانه‌ای هستند اقلیم 
آن‌ها وابستگی زیادی به میزان دما و بارش دارد و درنتیجه اين مناطق در برابر تغییر اقلیم دچار استرس 
می‌شوند (هیات بین‌الدول تغییر اقلیم؛ ۲۰۰۷). 

مدل‌های ]0/۷ به‌طور وسیع برای ارزیابی تغییر اقلیم در یک مقیاس جهانی استفاده می‌شود؛ اما خروجی 
این مدل‌ها برای ارزیابی تغییرات اقلیمی در سطح محلی و منطقه‌ای کافی و دقیق نیست. اکشر مطالعات 
ارزیابی تغییر اقلیم معمولا با دقت مکانی بالا روی مناطق کوچک انجام می‌شود. در این حالت شایسته است 
که توسط مدل‌های اقلیمی منطقه‌ای (186۷19) تا حد چندین ده کیلومتر یا کمتر ریزمقیاس گردند؛ اما ارتباط 
دادن تفکیک‌پذیری مدل‌های اقلیمی جهانی با مقیاس محلی یک فرآیند ریزاقلیمی است که خحود مسئله‌ای 
قابل توجه است. اخیراً تلاش‌های زیادی از طرف مجامع اقلیم‌شناسی برای توسعه روش‌های ریزمقیاس نمایی 
دینامیکی و آماری برای بیان تغییر اقلیم در یک مقیاس محلی و منطقه‌ای صورت گرفته است. مدل‌های 
71 نمونه‌ای از مدل ریزمقیاس نمایی دینامیکی است. مدل‌های دینامیکی منطقه‌ای دارای دقت فضایی 
مناسبی برای ارزیابی‌های اقلیمی بوده و همچنین فرآیندهای ریزمقیاس جوی و اقلیمی تا حد ممکن در آن‌ها 
لحاظ می‌شوند؛ اما مدل ]161۷ هنوز دارای چند اشکال است. برای مثال این مدل‌ها به مساحت مورد مطالعه 


و میزان تفکیک‌پذیری وابسته هستند و از لحاظ محاسباتی پیچیده بوده و انعطاف لازم را ندارند. مسئله 


۱۹۵ 
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دیگری که وجود دارد اینست که آموزش این مدل‌ها برای کاربران مدل و همچنین واسنجی آن‌ها بسیار 
هشتان یتشم ۲۱۳۳۱۵ 

به عنوان جانشینی برای مدل‌های منطقه‌ای دینامیکی. مدل‌های ریزمقیاس نمایی آماری توسعه بافته‌اند. این 
مدل‌ها بر اساس وضعیت اقلیم منطقه‌ای که نتیجه وضعیت اقلیم در مقیاس بزرگ است. پایه‌ریزی شده است؛ 
مانند توپوگرافی؛ توزیع دریا- خشکی و کاربری زمین (ویلبی ۲۰۰8). روش‌های ریزمقیاس نمایی آماری 
به‌طورکلی به سه طبقه تقسیم می‌شود: ۱. روش رگرسیونی؛ ۲. خوشه‌بندی الگوهای آب‌وهوایی و ۳. 
مولدهای سری‌های زمانی آب‌وهوایی. یکی از فواید بزرگ این روش‌ها در مقایسه با مدل‌های منطقه‌ای 
دینامیکی. ساده و ارزان بودن آن‌ها از لحاظ محاسباتی و کاربرد آن‌ها برای خروجی‌های مختلف مدل 007۷ 
برای سناریوها مختلف است. همین امر به ما امکان می‌دهد عدم قطعیت‌های سناریوهای آینده را ارزیابی 
کنیم. به‌طور نظری ضعف بزرگ ریزمقیاس نمایی آماری اینست که فرضیات پایه‌ای آن‌ها قابل تأیید نیست. 
یعنی روابط توسعه يافته برای اقلیم امروز برای شرایط اقلیم آینده پابرجا نیست (فولر۲۰۰۷). روش‌های 
ریزمقیاس نمایی آماری زیادی تاکنون توسط محققین هواشناسی و اقلیم‌شناسی ارائه شده است. دو روش 
استاندارد شده که کاربرد زیادی در ریزمقیاس نمایی مدل‌های تغییر اقلیم دارد و به‌راحتی توسط کاربران قاببل 
استفاده است. مدل 1۸1۹-۷۷6۶ که یک مدل تصادفی مولد سری‌های زمانی آب‌وهوایی و مدل ]915۷ که 
ترکیبی از مولد تصادفی آب‌وهوا و روش‌های رگرسیونی است. 

ویلبی و همکاران* (۲۰۰7) مدل‌سازی جامع اثر تغیبر اقلیم بر منابع آب و کیفیت آن برای رودخانه کنت 
انگلستان را ارائه نمودند. در این تحقیق, نتایج نشان داد که عدم قطعیت وسیع در سناریو تغییر اقلیم و منابع 
آب به علت انتخاب 001۷1 می‌باشد و این عدم قطعیت در فصل تابستان بسیار بیشتر است. 

سجاد خان و همکاران" (۲۰۰7) سه روش ریزمقیاس نمایی که عبارتند از مدل ریزمقیاس نمایی آماری 
(815[71). مدل 1۵135-۷۷0 و مدل شبکه عصبی مصنوعی (((۸) را از طریق آزمون‌های عدم قطعیت 
برای سه پارامتر بارش روزانه. حداقل دمای روزانه و حداکثر دمای روزانه مورد مقایسه قرار دادند. نتایج 
نشان داد که روش ٩251۷‏ بهترین کارایی و روش ۸ کمترین کارایی را دارد و روش 1۵15-۷70 در 


حد وسط قرار دارد. 
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کارآموز و همکاران" (۲۰۰۹) روش ریزمقیاس نمایی ۹1251۷ و روش شبکه‌های عصبی ((۸۵) را برای 
پیش‌بینی بلندمدت بارش در ۵ ماه از سال (دسامبر تا آوریل) در جنوب شرقی ايران بکار بردند. نتایج نشان 
داد که روش ]8125[۷ عملکرد بهتری نسبت به روش ۸ دارد. ساروار و همکاران" (۲۰۱۰) اثرات تغییر 
اقلیم را بر حوض آبریز رود تایمز با استفاده از مدل‌های ریزمقیاس نمایی ۹15 1۵16-۷۷0و 
( 6 مورد ارزیابی قرار دادند. نتیجه اينکه مدل 9129/۷ برای دماهای حداقل و حداکثر عملکرد 
بهتری نسبت به دو مدل دیگر دارد ولی برای بارش مدل 1185-۷۷ ترجیح داده می‌شود. 

مینو و همکاران" (۲۰۱۲) از یک مدل هیدرولوژیکی (3.4 1۳0-111۷5 برای مدل بندی هیدرولوژیکی 
منطقه رودخانه تونگا-بهادرا در هندوستان و از روش ریزمقیاس نمایی 1251۷ برای ریزمقیاس نمایی 
دماهای حداقل و حداکثر و بارش روزانه در منطقه مورد مطالعه استفاده کردند. نتایج مطالعه بیلان آب؛ 
افزایش بارش و رواناب و کاهش تلفات تبخیرتعرق واقعی را روی منطقه مورد مطالعه نشان داد. 

رجبی و شعبانلو* (۲۰۱۲) از مدل 1 برای ارزیابی تغییر اقلیم منطقه کرمانشاه در غرب ایران و اثر 
آن بر شاخص‌های اقلیمی همچون جانسن. کرنر و دومارتن استفاده کردند. نتایج نشان داد که اقلیم منطقه 
بخصوص در دورهٌ ۲۰۷۰-۲۰۹۹ خشک‌تر خواهد شد و اين تغییر برای سناریوی ۸2 محسوس‌تر است. 

چیما و همکاران* (۲۰۱۳) عملکرد روش ریزمقیاس نمایی 5151۷ را روی روند داده‌های حداقل دمای 
ایستگاه‌های کشور پاکستان را در دورة ۱۹۹۱-۲۰۱۰ مورد ارزیابی قرار دادند. نتایج نشان داد که طبق آزمون 
من کندال روند افزایشی حداقل دمای سالانه معنی دار است. نتایج نشان داد که تغییر معنی دار اقلیمی بیشتر 
در نواحی شمالی کشور پاکستان رخ داده است. کاظمی و همکاران" (۲۰۱۶) روش ٩1251۷‏ را روی خروجی 
دمای روزانه مدل جهانی ۳1۸۱۷5 بکار بردند. نتایج نشان داد که داده‌های ریزمقیاس شده خیلی دقیق‌تر 
از داده‌های مدل ۲۳۱۸۷5 است. 

نوری و آلام" (۲۰۱۶) مدل ٩1251۷1‏ را روی داده‌های مشاهداتی دما و بارش روزانه از سال ۱۹۸۱ تا 
با استفاده از مدل جهانی ۳1۸00۷13 بکار برد و نتایج نشان می‌دهد که داده‌های دما و بارش ریزمقیاس 


شده با داده‌های مشاهده شده تقریباً مطابقت داشته‌اند. 


1 12۲210007 20 6۲ 1 

لچ 1 20 92۳7۷۵۲ 2 

لح 6 20 ۱۷6۵80 3 
۲امصحطاقدی 0صم .1طوزف۴ 4 
2 ۲ 200 قصععط) 5 

6 122101 20 61 1 

مداخ مضه .1 ۸ لا 7 


سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی... ۷ 


محمدی و همکاران (۱۳۸۹) اثر تغییر اقلیم را بر پارامترهای هواشناسی و آبشناحتی حوضه آبریز قزل 
اوزن در استان زنجان بررسی نمودند. آن‌ها رژیم ماهانه بارش و دما حاصل از مدل ۲1۸001۷13 تحت دو 
سناریوی اقلیمی ۸2 و 81 را با داده‌های مشاهداتی مقایسه کردند. نتایج پیش‌بینی تقریبی پارامترهای 
هواشناسی و هیدرولوژیک تا سال ۲۰۵۷ میلادی را به صورت کاهش میانگین بارش سالانه به میزان حدود ۲۱ 
میلیمتر و افزایش میانگین دمای سالانه حدود ۲/۵ درجه سانتیگراد و کاهش میانگین سالانه رواناب حدود ۱۰ 
متر مکعب در ثانیه نسبت به میانگین دوره مشاهداتی (۱۹۷۲-۱۹۹۹) نشان داد. 

دهقانی پور و همکاران (۱۳۹۰) در مقاله خود توانمندی مدل ریزمقیاس نمایی آماری ]81«5[۷ را در تولید 
داده‌های بارش, دما و تبخیر در ایستگاه هواشناسی سینوپتیک تبریز مورد بررسی قرار دادند. نتایج حاصل 
نشان می‌دهد که مدل 1291۷ توانایی مناسبی در کوچک مقیاس نمودن داده‌های دماء تبخیر و بارش دارد. 

گل محمدی و مساح بوانی (۱۳۹۰) در مقاله خود به ارزیابی اثرات تغییر اقلیم بر وضعیت خشکسالی 
حوضه قره سو در دوره‌های آتی با استفاده از شاحص خشکسالی بارش استاندارد شده ([8۳) پرداختند. در 
این ارزیابی از روش ۹1251۷۲ برای ریزمقیاس نمایی خروجی مدل جهانی 112007۷3 تحت سناریوی ۸2 
استفاده شده است. مقایسه نتایج حاکی از این است که حوضه در دوره ۲۰۶۰-۲۰3۹ میلادی شاهد افزایش 
بارش و به تبع آن کاهش شدت خشکسالی نسبت به دوره پایه خواهد بود. این کاهش شدت. برای مقیاس- 
های زمانی + و ۱۲ ماهه با افزايش دوره بازگشت افزایش می‌یابد. این در حالی است که برای خشکسالی های 
۶ ماهه با افزايش دوره بازگشت شدت خشکسالی حوضه بیش از دوره پایه خواهد شد. 

با توجه به مطالعات فوق مشاهده می‌شود که مدل ٩1951۷‏ دارای دقت قابل قبولی در ریز مقیاس نمایی 
داده‌های اقلیمی می‌باشد گرچه بیشتر مطالعات صورت گرفته بر اساس سناریوهای" (۸4) می‌باشد ولی 
تحقیق حاضر در نظر دارد که با استفاده از حروجی 2 که یکی از مدل‌های تغییر اقلیم جفت 
شده "01۷11۳5 و بر اساس سناریوهای جدید" (۸۵125) می‌باشد. متغیرهای اقلیمی دما و بارش را در طول 
قرن ۲۱ شبیه‌سازی و مورد بررسی قرار دهد. 
۲ منطقه مورد مطالعه 

تبریز مرکز استان آذربایجان شرقی با مساحتی حدود ۱۲۰۰ کیلومتر مربع در ۶7 درجه و ۱۷ دقیقه طول 


شرقی و ۳۸ درجه و ۰۵ دقيقه عرض شمالی واقع شده است. ارتفاع آن از سطح دریا ۱۳۹۲ متر است. راکو 


۱ گزارش چهارم 1۳16 
۲ فاز ۵ پروژه متقابل مدل جفت شده 


۳ گزارش پنجم)1۳ 


۱5۸ جغرافیا و مخاطرات محبطی شماره پیست و یکم 


از شمال به کوه عینالی و از جنوب به دامنه کوه سهند و از سمت غرب به جلگه تبریز و دریاچه 
ارومیه محدود است. شکل (۱) موقعیت محدوده مورد مطالعه را در سطح استان و کشور نشان می‌دهد. 


شکل ۱ موقعیت جغرافیایی ایستگاه تبریز در کشور و استان آذربایجان شرقی (نگارندگان) 


۳ مواد و روش‌ها 
در این تحقیق برای بررسی اثر تغییر اقلیم در محدوده مورد مطالعه از ایستگاه تبریز به دلیل داشتن آمار 
بلند مدت استفاده شد. دوره آماری مشاهداتی پایه دوره ۳۰ ساله ۱۹۱۱-۱۹۹۰ است که در جدول شماره (۱) 
آورده شده است. داده‌های روزانه حداقل دماء حداکثر دما بارش در دوره مربوطه وارد مدل شده و پس از 
ارزیابی کیفی و آماری توسط نرم افزار و همچنین تأیید اعتبار مدل از طریق انجام آزمون‌های آماری داده‌ها 
جهت تولید سری‌های تصادفی روزانه آماده شدند. 
جدول ۱ اطلاعات مربوط به ایستگاه سینوپتیک تبریز 


ارتفاع از سطح دریا عرض جغرافیایی طول جغرافیایی طول دوره آماری 


(متر) (درحه) (درحه) پایه 


ردیف | ایستگاه 


۱ تیار ۱۳۹۱ ۱۳۸۹ 1-۷ ۱۹۱-۰ 
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در این تحقیق از حروجی مدل 20۳5[۷2» که ازچهارمین نسل از مدل‌های آب‌وهوایی است که توسط 
مرکز مدلسازی و تحلیل آب‌وهوای کانادا (60070) زیر نظر سازمان محیط زیست این کشور توسعه یافته 
است. در این مدل کل سطح زمین بصورت ۱۳۸۶ سلول. شبکه بندی شده است که مشخصات آن در 
جدول ۲ ارائه شده است (6021708,2014). 


حدول ۲ مشخصات مدل 6206591۳02 


دقت مکانی داده‌های اقیانوسی دقت مکانی داده‌های اتمسفری 
طول جغرافیایی عرض جغرافیایی طول جغرافیایی عرض جغرافیایی ۳ 
(درحه) (درحه) (درحه) (درحه) 
۱۰۵ ۷ و ۰/۹۳۰۳ ۲۱۳۸۵ ۳۸/۹۰۹ 072 


۰ ۳. سناریوهای میزان غلظت گازهای گلخانه‌ای 

هیات بین الدول تغییر اقلیم در تدوین گزارش پنجم ارزیابی خود (۸15) از سناریوهای جدید 1۳ به 
عنوان نماینده‌های خحطوط سیر غلظت‌های گوناگون گازهای گلخانه‌ای استفاده کرده است. سناریوی جدید 
انتشار دارای چهار خحط سیر کلیدی با نام‌های130۳2.6, 130۳45 ۳6)لو 6۳8.5 می‌باشند که بر 
اساس میزان واداشت تابشی آن‌ها در سال ۲۱۰۰ نام گذاری شده‌اند (627700,2014). 

- سناریوی انتشار 13۳2.6 

این سناریو توسط تیم مدل‌سازی 11۸01 از موسسه ارزیابی‌های محیطی هلند طراحی شده است. این 
سناریو در برگیرنده کمترین نرخ افزایش گازهای گلخانه‌ای و واداشت تابشی ناشی از آن است. مطابق این 
سناریو واداشت تابشی در اواسط این قرن به حدود ۳,۱ وات بر متر مربع رسیده و سپس کاهش يافته و به 
وات بر مترمربع در سال ۲۱۰۰ می‌رسد. برای رسیدن به این سطح واداشت تابشی بایستی گازهای 
گلخانه‌ای به میزان قابل توجی کاهش يابند. 

- سناریوی 1304.5 

سناریوی 16۳4.5 توسط گروه مدل‌سازی 1۷11010۵۷۲ طراحی شده است و در آن واداشت تابشی 
ناشی از گازهای گلخانه‌ای قبل از سال ۲۱۰۰ در مقدار ۶,۵ وات بر مترمربع ثابت می‌ماند. 

- سناریوی 13)۳6 

سناریوی انتشار 136۳6 توسط گروه مدل‌سازی ]۸1۷ در موسسه ملی مطالعات محیطی ژاپن طراحی 
گردید. در این سناریو واداشت تابشی بعد از سال ۲۱۰۰ به دلیل استفاده از فناوری‌های جدید و سیاست‌های 


کاهش گازهای گلخانه‌ای ثابت می‌ماند . 
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- سناریوی 10۳8.5 

بدون اتخاذ هیچ گونه سیاست‌های کاهش آثار و مقابله با پیامدهای اقلیم. آب‌وهوای کره زمین در حط سیر 
سناریوی انتشار 8.5 )1 پیش خواهد رفت. به طوریکه ادامه این روند منجر به واداشت تابشی به میزان ۸,۵ 
وات بر متر مربع در سال ۲۱۰۰ می‌گردد. در این هنگام غلظت دی اکسید کربن به ۱۰۰۰ 0۳10 رسیده و همچنان 
روند افزايش خواهد داشت. این سناریو توسط تیم مدل‌سازی ۷۳۹5۸1۲ و موسسه 11۸5۸ به سرپرستی 
پروفسور کیوان ریاحی در موسسه بین المللی آنالیز سیستم‌های کاربردی 11۸5۸ اتریش توسعه و طراحی شد که 


وجه مشخصه آن روند افزایش گازهای گلخانه‌ای است. 


۲ ۳ ریز مقیاس کردن داده‌های اقلیمی 

در این مطالعه از روش آماری جهت ریز مقیاس کردن خروجی مدل‌های 06۷1 استفاده می‌شود. تاکنون دو 
رویکرد برای روشهای ریزمقیاس نمایی آماری ارائه شده است. رویکرد اول که بیشتر مورد استفاده قرار می‌گیرد 
استفاده از مدل‌های رگرسیونی یا تابع انتقال است که در آن با استفاده از رابطه بین پارامترهای مختلف جوی و 
سینوپتیکی (متغیرهای پیش‌بینی کننده) و پارامتر اقلیمی که لازم است دورنمایی از آن در آینده داشته باشیم (متغیر 
پیش‌بینی شونده) یک تابع انتقال ارائه می‌شود. رویکرد دوم استفاده از مولدهای تصادفی آب وهوایی است که بر 
اساس مدل‌های سری زمانی» زنجیره‌های مارکف و سری‌های فوریه ارائه می‌شود. یکی از روشهایی که بر اساس 
این رویکرد بنا شده است. مدل آماری 1۵6-۷۷0 است (1991 ,له 66 16910 یکی از نرم افزارهایی که 
براساس ترکیبی از این دو رویکرد پایه ریزی شده ٩1951۷۲‏ است. این نرم افزار داده‌های روزانه هواشناسی 
ایستگاه مورد نظر را به عنوان ورودی گرفته و در هفت مرحله خروجی مدل ]20۷) را بر اساس داده‌های روزانه 
در منطقه مورد نظر ریزمقیاس می‌نماید. 

این هفت مرحله عبارتند از: 

۱. کنترل کیفی و تبدیل داده‌ها: ایستگاه‌های هواشناسی کمی هستند که دارای داده‌های بدون نقص باشند. نذا 
شناسایی داده‌های غلط یا گم شده ضروری است. این نرم افزار یک چک کنترل کیفی ساده در ابتدا روی داده‌های 
روزانه ایستگاه انجام می‌دهد. علاوه بر آن گاهی اوقات لازم است داده‌ها جهت برازش بهتر مدل تبدیل شوند 
(برای مثال تبدیل لگاریتمی توانی» وارونگی و غیره). 

۲.بررسی و انتخاب متغیرهای پیش‌بینی کننده: از آنجا که اساس کار این روش تعیین رابطه بین متغیرهای 
پیش‌بینی کننده سینوپتیکی (ازقبیل فشار سطح دریا) و متغیر پیش‌بینی شونده يا پارامتر اقلیمی منطقه (از قبیل بارش 


ایستگاه) است در این مرحله متغیر پیش‌بینی کننده‌ای که دارای همبستگی بیشتربا متغیر پیش‌بینی شونده دارد و یا 
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سهم بیشتری در تغییرات مدل دارد بررسی و انتخاب می‌شود. البته در این مرحله تجربه و دانش کارشناس 
هواشناسی اهمّیت زیادی دارد. 

۳ برازش مدل: در اين مرحله مدل رگرسیون چند متغیره براساس متغیرهای پیش‌بینی کننده انتخابی در مرحله 
قبل و متغیر پیش‌بینی شونده به روش حداقل مربعات ساخته می‌شود. کاربر ساختار مدل را بر این اساس که 
ماهانه» فصلی و پا سالانه باشد. تعیین می‌کند. علاوه بر آن باید مشخص شود که مدل شرطی است يا غیر شرطی. 
در مدل غیرشرطی فرض می‌شود که متغیر پیش‌بینی شونده مستقیماً با متغیرهای پیش‌بینی کننده ارتباط دارد. ولی 
در مدل شرطی یک فرآیند بینا بين بین عوامل منطقه‌ای و هوای محلی وجود دارد (برای مثال مقدار بارش محلی 
به رخدادهای خشک و تر وابسته است که به‌طور چرخشی وابسته به پیش‌بینی کننده‌های مقیاس محلی از قبیل 
رطوبت و فشار هواست». 

؛.تولید داده‌های آب‌وهوایی: مولدهای هواشناسی نمای کلی از سری داده‌های مدل شده روزانه بر اساس 
متغیرهای پیش‌بینی کننده مشاهده شده (یا داده‌های درون یابی شده ۲ تولید می‌کند. این پروسه قادر به 
تأیید مدل برازش شده (با استفاده از داده‌های مستقل) و سری‌های زمانی مدل شده برای شرایط اقلیمی حاضر 
است. در این مرحله کاربر می‌تواند طبق هر دوره زمانی داده هواشناسی تولید کند و آن‌ها در مراحل بعد با 
داده‌های واقعی مورد مقایسه قرار داده و اعتبارسنجی مدل را انجام دهد. 

۵تحلیل داده‌ها: در این مرحله می‌توان میانگین‌ها. تحلیل فراوانی وشاخص های آماری داده‌های مشاهده شده 
و سناریوهای تولید شده را محاسبه و تحلیل قرار داد. در این بخش کاربر باید دوره زمانی مربوطه و نام فایل 
خروجی را مشخص نماید. 

1.تحلیل گرافیکی: در این مرحله تحلیل فراوانی داده‌های به روش‌های نموداری انجام می‌شود که در آن توزیع 
اماری داده‌های تولید شده و مشاهده شده مورد تحلیل قرار می‌گیرد و داده‌های پرت و حدی مشخص می‌شود. 
این تحلیل‌ها شامل توزیع‌های آماری تجربی, گامبل. نمایی و حدی است. علاوه بران در این مرحله میانگین‌های 
ماهانه وفصلی داده‌های دوره‌های مختلف با هم مقایسه می‌شوند. 

۷تولید سناریو: در نهایت عملیات تولید سناریو بر اساس داده‌های روزانه تولید شده مدل انجام می‌شود. در ایین 
مرحله خروجی مدل جهانی اقلیمی را در نظر گرفته و بر اساس داده تولید شده مدل در دوره‌های آینده خروجی مدل 


ور انتخاب متغیرهای پیش پینی کننده 
تغییرات و ویژگی‌های پارامترهای ایستگاهی مورد نظر را توصیف کند؛ که این متغیرها باید از متغیرهای )۷ 
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انتخاب شود. متغیرهای ۷2 شامل ۲۱ متغیر اتمسفری است و از میان آن‌ها متغیر یا متغیرهایی انتخاب می- 
شوند که دارای بالاترین میزان همبستگی و کمترین مقدار شاخص *1 می‌باشد و به عنوان متغیرهای توضیحی 
جهت پیش‌بینی پارامترهای دما و بارش وارد معادله رگرسیونی مدل ٩1251۷1‏ می‌شوند. جهت محاسبه ضریب 
همبستگی از مدل ٩1251۷]‏ و جهت محاسبه شاخص ۳ از رابطه زیر استفاده شده است: 

2, 

حت ۳ ۱ 


۶ ۳ ارزیابی و واسنجی مدل 
به منظور کالیبره نمودن و واسنجی مدل ٩1251۷]‏ داده‌های مورد استفاده, ایستگاه مشاهداتی شهر تبریز و 
داده‌های مرکز ملی پیش بینی متغیرهای محیطی کانادا ۲0۳۳ به دو دور 15 ساله (۱۹۲۰-۱۹۷۵) و (1۹9۰- 
۲ تقسیم شدند که از 15 سال اول برای کالیبره کردن مدل با استفاده از روش بهینه سازی حداقل مربعات 
استفاده شد. بعد از کالیبره نمودن مدل به منظور حصول اطمینان از این که مدل توانایی شبیه سازی داده‌ها را خارج 
از محدودهة زمانی کالیبراسیون دارد. نیاز است مدل ارزیابی شود بدین منظور با استفاده از مدل ٩1۸9/۷]‏ کالیبره شده 
دما و بارش را برای دوره (۱۹۰-۱۹۹۰) شبیه سازی شدند. سپس با مقايسة داده‌های مشاهده‌ای و داده‌های شبیه 
سازی شده کارایی مدل برای ایستگاه تبریز مورد بررسی قرار گرفت .در این تحقیق به منظور بررسی میزان کارایی 
مدل از دو معیار آماری میانگین حطای مطلق (۱)0۷۲۸7. درصد اریبی (۲)۳۳1۵5 استفاده گردید. این معیارها بر 
اساس روابط ۲ و ۳ محاسبه می‌شوند. 
0( - :۱0۵ 
)سر 


که در آن: 0 داده‌های مشاهداتی» ,5 داده‌های برآورده شده می‌باشد. مقادیر ۱/۳۹۹ معرف دقفت در بسیاری 


)۲( ۸/۸ 21/۸۷2 )6,-0( )۲( ۳8/۸۹ < 


از روش آماری می‌باشد. هر چه این مقادیر به صفر نزدیک‌تر باشد نشان دهنده دقت بالای مدل و مقدار صف 
نشان دهنده عدم وجود خطا در برآورد مدل است. ۳191517 هرچه به صفر نزدیک‌تر باشد نشاندهنده دقت بالای 
مدل در برآورد متغیر مورد نظر می‌باشد. در صورتی که مقدار این پارامتر به سمت مثبت میل کند نشان دهنده ایین 
است که متغیر مورد بررسی بسیار کمتر از مقدار واقعی برآورد شده است و بر عکس. قابل ذکر است که برای این 
پارامتر آستانه خاصی در نظر گرفته نشده است (211,2007 66 0۷1011251 

در انتخاب پیش‌بینی کننده‌ها نرم افزار 51081۷1 در واقع بین سری مشاهداتی روزانه منطقه و متغیرهای بزرگ 
مقیاس مشاهداتی منطقه (۳۳() رابطه‌ای با حداکثر ضریب همبستگی برقرار نموده (گل محمدی و مساح 


۱۱۱ 
عقلط 06۲۵۵۳۲ 2 
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بوانی۱۳۸۹۰) و سپس با استفاده از پارامترهای بدست آمده از این رابطه با بکارگیری متغیرهای بزرگ مقیاس 
حاصل از مدل 0210۳5[۷12 برای دوره ۱۹۲۰-۱۹۹۰ استخراج شدند. سپس با معرفی این داده‌ها به مدل 51951۷1 
بارش و دمای حداکثر و حداقل تحت سه سناریو 1002.6,۲004.5,7008.5 برای سه دوره ۳۰ ساله ۲۰۱۰-۲۰۳۹ 
و ۲۰۰-۲۰7۹ و ۲۰۷۰-۲۰۹۹ برای ایستگاه مورد نظر تولید کرده و تغییرات مقادیر شبیه‌سازی شده در 


سناریوهای مختلف نسبت به دوره پایه مورد بررسی قرار گرفت. 


۶ بحث و نتایج 


۶ انتخاب متغیرهای مستقل جهت شبیه‌سازی و کالیبره کردن مدل 

در مدلهای تابع انتقالی مانند مدل ٩1951۷1‏ قبل از کالیبره کردن مدل لازم است متغیرهای اقلیمی که بیشترین 
همبستگی را با پارامترهای مورد نظر دارند. تعیین گردند. در پژوهش حاضر ابتدا با استفاده از پارامترهای بزرگ 
مقیاس مشاهداتی ۱۱0۳۲ و برم افزار 71 به انتخاب پیش‌بینی کننده‌های منتخب پارامترهای اقلیمی مورد نباز 


اقلیمی مورد نظر انتخاب گردید که نتایج پیش‌بینی کننده‌های منتخب نهایی در جدول ۳ ارائه شده است. 


حدول۳ انتخاب متغیرهای پیش‌بینی کننده بزرگ مقیاس برای پیش‌بینی پارامترهای محلی بارش. حداقل دما و 
حداکثر دما در ایستگاه تبریز 


۲ نگ ۳ 5 
متغیر پیش‌بینی کننده‌های منتخب ۱ ۱ + | ۳۲۵۱۲۵ 
مطلق (ج1) | جزئی (7) 
ميانگین فشار سطح دریا ۹ ۹- ۳۳ ۰/۰ 
حداقل رطوبت نسبی در ارتفاع ۰ هکتوپاسکال ۰۸۳۹ ۰/۳۹ 2 ۱/۰ 
دما رطوبت نسبی در ارتفاع ۵۰۰ هکتوپاسکال 2 ۰۸/۸۹ ۱/۱ ۱/۰ 
میانگین دما در ارتفاع ۲ متری ۰/۷ ۳۵ اأ ا ۵ ۰/۰ 
تاوایی در ارتفاع ۵۰۰ هکتوپاسکال ۰/۱۸۰ 2۳ ۰/۸ ۰/۰ 
حداکثر رطوبت نسبی در ارتفاع ۸۵۰ هکتوپاسکال ۸۸-۳ ۸۸۲ ۰/۰ ۰/۰ 
دما رطوبت نسبی در ارتفاع ۰ هکتویاسکال ۰,۳۲ ۰/۸۱۸ ۱/۸ ۰/۰ 
ميانگین دما در ارتفاع ۲ متری ۷۰۵ ۰/۹۹ ۳3۹ ۰/۰ 
ملفه سرعت نصف النهاری در سطح ۵۰۰ هکتوپاسکال ۰/۳۵ ۰/۱۱۳۱ :1/1 ۰/۰ 
تاوایی در ارتفاع ۵۰۰ هکتوپاسکال ۷ ۰۹۶ ۰/۶۰ ۰/۰ 
بارش واگرایی در ارتفاع ۵۰۰ هکتوپاسکال ۰:۹ ۰۱۵ ۲ 2 
واگرایی در ارتفاع ۰ هکتوپاسکال ۰۱۳۲ 2 ۰/۹ ۱/۰۰۵ 
رطوبت نسبی در ارتفاع ۵۰۰ هکتوپاسکال ۱:۷ ۸۳ :1 ۰/۰ 


۵۵۵ ۳۱۷1۳0۳۱۵۵ ماع ۱۲2010821 1 


۱۳۹۶ جغرافیا و مخاطرات محبطی شماره پیست و یکم 


با توجه به نتایج جدول ۲ مشخص شد که داده‌های دما هنت کر بهتری با داده‌های مشاهداتی داشته‌اند این 
بدان علت است که تغییر پذیری دما نسبت به بارندگی کمتر است و دما یک متغییر پیوسته می‌باشد که کمتر 
تحت تأثیر ناهنجارهای زمانی قرار می‌گیرد. در حالی که بارش تحت تأثیر عوامل مختلفی قرار می‌گیرد. از سوی 
دیکره بارند کی یک پارامتز کسشته انستا: 


۲ ۶ نتایج حاصل از ریز مقیاس نمایی توسط رم افزار ۱( 
در این بخش با بررسی و مقایسه داده‌های مشاهداتی و داده‌های شبیه‌سازی شده دما و بارندگی 7 ها 
در دوره پایه (۱۹۳۲۱-۱۹۵۹۰) برازش مدل بر اساس متغیرهای پیش‌بینی کننده ۱65 مورد بررسی و ارزیابی 


قرار گرفت و برای روشن شیلان توانمندی مدل علاوه بر روش ترسیمی از معیارهای آماری ( ۷۲۵,۴۳۸ 
نیز استفاده شد. نتایم حاصل از ارزیابی مدل با استفاده از معیارهای آماری در جدول ۶ ارائه شده است. 


جدول ۶ معیارهای آماری ارزیابی کارایی شبیه‌های ]1251۷ ایستگاه تبریز طی دوره آماری پایه (۱۹۰-۱۹۹۰) 


حداقل دما حداکثر دما بارش 

شاخص ۳۳۵۹ ۱/۳۰۹ ۳۳4۵۵ ۱/۳۰۹ ۳۵ ۱/۳۰۹ 
ژانویه ۷- | ۰,۲۹۲۳۵۸ ۷ - | ۱۱۲۳۵ [ ۲,۰۵۸ 
فوریه ۷ ۱ ۳ [ ۱۱۲ | ۳۷۵ و۱۵۵ ۳۸۷۶۹۵ 
مارس ۱۱۲( ۱ - 9۹۹۹« ۱« ۹ ۱۳ ۱ 
آوریل ۳ ۲ - ۰,۸ - ۳ ۰۱ ۱9۹9۳۷۵ 
می ۶ - ۵ ,- ۱:2۷ "ص«,۰ ۱۲ ,۰ ۹ ۸ ۵ ۲( 
ژوئن ۳۵ ۰,9۳۸ اک | ۰,۱۹۵ | ۱۵,۷۷۳۲- ۳۳۹ - 
جولای ۸,- ۵۹ | ۰۱۱۹۲۳ ۰۰۳۱۹۶« ۱ ۱ 
آگوست ۱۷۸ ۱۹-۸ ۳ | ۲۳ [ ۵۱۰۷۳ - ۱ 
سپتأمبر ۱/1 ۰,9۱۸ 99۰ص« ۰,۱9۹,۱۳۸ ۹,- ۷۹,- 
اکتبر ٩1۷‏ ۰/۰۹۳ ۱ ۳۵ ۱۲۳ ۱ 
توامیر ۹۸ ۱۳۳ ۱ص« ۳ ۱ - 1+۷" ۱۳ 
دسامبر ۲ - 9۹ ۸ | ۰,۰۱۵۱۸- ۱۳ ۱۳ 
زمستان ۳ ۸۶" ۱- ۳ و ۳ - ۷ 
بهار ۹,- ۲ - ۹۹ ۱ و۱۳ و 
تابستان ۲ ۲ ۳۱ ۰۲ ۱ - ۳۵ - 
پائیز ۱,۶ ۰:9۵ ۲۱ - ۲ - ۹ ۳ 
سالائه ۰,۱۶۹ ۸۲ | ۰,۰۳۸ 1 ۲,۱۳۸۷۹- ۱ 


سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی... 


دما (درجه سانتی گراد) 
دما (درجه سانتی گراد) 


1[ یج 
۵ تّ ف ند ئ 
بویتوی 
ی کي تسه کي 
که گی 
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شبیه سازی شده 7 77*7 مشاهداتی 7*۳ 


شکل ۳۲ داده‌های مشاهده شده و داده‌های مدل شده شکل۳ داده‌های مشاهده شده و داده‌های مدل شده 


میانگین دمای حداقل طی دوره (1۹7۰-۱۹۹۰) میانگین دمای حداکثر طی دوره (۱۹۰-۱۹۹۰) 


بارش (میلیمتر) 


کي کی ی خی ٍ و «ج ند ی دی 
کی کی کی کی ات کب کلب کم کي کی کي 
کی کتی کي کم کي اي 
هي ۰ 


شبیه سازی شده -ف مشاهداتی سم 


شکل ۶ داده‌های مشاهده شده و داده‌های مدل شده میانگین بارش طی دوره (۱۹7۰-۱۹۹۰) 


با توجه به نتایج جدول؛ و شکل‌های ۲ تا ؛ مشخص شد که داده‌های شبیه‌سازی شده توسط متغیرهای 
پیش‌بینی کننده منتخب با داده‌های مشاهداتی در طول دوره پایه (۱۹۱۰-۱۹۹۰) تطابق بیشتری دارند که در ایین 
میان تطابق دمای حداکثر مشاهده شده با شبیه‌سازی شده بیشتر از پارامترهای بارش و دمای حداقل می‌باشد؛ 
بنابراین بر این اساس مشخص کگردید که متغیرهای پیش‌بینی کننده منتخب از کارایی مناسب و قابل قبولی 
جهت شبیه‌سازی پارامترهای اقلیمی مورد نظر برای دوره آتی دارند. 
۳ . شبیه سازی پارامترهای اقلیمی ایستگاه تبریز تا سال۲۱۰۰ 

در این بخش پس از ارزیابی دقت مدل برای دوره پایه (۱۹۱۱-۱۹۹۰) اقدام به شبیه‌سازی پارامترهای 
اقلیمی دما و بارش برای دوره‌های آتی شد. مدل را برای سه دوره ۳۰ ساله در قرن ۲۱ (۲۰۱۰-۲۰۳۹ ۲۰۹- 
۰ و ۲۰۷۰-۲۰۹۹) و براساس مدل جهانی 620۳91۷12 تحت سه سناریو 1002.6,1004.5,۲0۳08.5 برای 


ایستگاه تبریز تولید و تغییرات مقادیر شبیه‌سازی شده در سناریوهای مختلف نسبت به دوره پایه مورد بررسی 


قرار گرفت. 


۱۳۹۹ جغرافیا و مخاطرات محبطی شماره پیست و یکم 


- بررسی تغییرات بارش 

جدوله تغییرات بارش را تحت سناریوهای 100 برای سه دوره نسبت به دوره پایه نشان می‌دهد. جهعت 
مقایسه. نمودار داده‌های مربوطه در شکل‌هایه آورده شده است. باتوجه به جدول مشاهده می‌کنیم که 
به‌طورکلی بارش در سه سناریوی مورد بررسی برای دو دوره ۲۰۱۰-۲۰۳۹و ۲۰۷۰-۲۰۹۹ کاهش بارش و 
برای دوره ۲۰۶۰-۲۰۹۹ افزايش بارش خواهد داشت. همچنین بارش به‌طورکلی در فصل زمستان افزایش و بقیه 
فصول با کاهش بارش مواجه خواهد بود. تغییر اقلیم به دو شکل می‌تواند بر روی بارش نقش داشته باشد. یکی 
به شکل تغییر در مقدار بارش و دیگری تغییرات زمانی الگوی بارش می‌باشد. نکته جالب توجه در این تحقیق 
تغییرات زمانی الگوی بارش می‌باشد که هرچه به پایان قرن حاضر نزدیک می‌شویم ملاحظه می‌گردد که بارش 
در ماه‌های اول سال کاهش پیدا کرده و در ماه‌های آحر سال, افزایش پیدا می‌کند. 


جدول ۵ تغییرات بارش تبریز طی دوره ۲۰۷۰-۲۰۹۹ و۲۰۱۰-۲۰۶۰,۲۰۳۹-۲۰۹۹ در مقایسه با دوره پایه ۱۹۲۰-۱۹۹۰ 


دوره 13۳26 1+05 رز ۱4۵ 
۶ گ ّ ک گ 4 ى 
ستاریو | زا 1 1 / ۲ / 1 ۲ 
2 ‌ 2 ۳۹ ی 5 
ء ح ب ح تء بِ حء حء بِ 
ژانویه ٩,۰۶۳۹۲-‏ -۲۱,۵۳۱۱ ۱۳۶ ۱۰,۱۵۹ -۶۲,۵۵۵۱ ۱۳۹۹۰۵ ۱۳۱۲۸۰۲ ۱۵,۶۳۷ 3۳3222 
فوریه ۳9 +۳ ۱ 33 ۱۳۸ ۹( 0۶۳ ۱۳۹4 ۸۵۰ 
مارس ۲ ۶۲,۱۱۸ ۱۶۸۷/۸۹ ۱۳/۶ -۳۲,۱۲۱۲ ۱۷۱۲۵۸۰۲ ۸۹ -۶۲,64۵ ۱۳۸ 
آوریل ۰۱۲۶۲ ۱ ۱۳۸ ٩2۷ ۱ ٩,۱۱۲‏ ۳۹۹۹ 1۱ ۱ ۱0۳ 
می ۲۱,۱۸۱ | ۵۵,۱۵۷ -۵۳۲ ,۶۸ ۳۹/۵۱ ۵۳۱۷۵۸۵ ت۱9 ۹ ۳۸۵ 
ژوئن ۸ -۱۹,۹۸۲۱ ۱۱۳۳ ۰ ۶۸۸۶۱۵ -۱۳,۹۹۱ ۵۹/۰۹۳ ۸۳۹۹۲ ۲۸۸۹۸۱ 16,۱۷۹ 
حولای ۶۰۳۸۵۲ 0۳ 1۲,۸۷۱ ۳۷,۲۸۵ ۱۹,۹۵۱ 1۳,۶۵۱ ۲۹,۰۱۹ ۰۲۷۶۷۸ 1۱۲۰۹۵۹ 
آگوست | -۵۸۳۲۲ | ۹۵,۳۵۰ ۷۰۳۷/۲۱ ۷,۱۵۲ 2 ۱۷۲۱۲۹۵ 1 -۸۵۸۲۲ ۱۸ 7۱۷۲۷۸۷ 
سیتامبر 14,۱7۳ ۱۷ ۲۶,۲۷۰ ۵۱۱۸۸۷۹ ۱ ۶۷۷۵۳۲ ۳۵,۵ ۱۷۸۹۶ ۵۸,۶۲۶ 
اکتیر ۳۹,۱۸۲ ۱۶,۹ 13-۳ ۶۱,۶۹1 ۱۱۷۶ ۹ ۲۸,۱۷۷ ۹,۶۱ 3۳۹4/۸ 
نوامبر ۱ ۸۳۹ ۱۷ ۹۵ ۱/۹ ۱,۷ ۱۳۹۳ ۲۲,۸۳۶ ۳ ۱ ۱20 
دسامپر ۸۵ ۷۲ ۹۳/۱ 2-9۰۶" ۱۲۱۳۹۳ ۳ ۱( ۱4 ۱ ۸۰۳ 
زمستان ۱۸۳ ۲ ۱ 3 ۱۳۱۳ ۸۹۶ 0۳۷ ۱۱۳ ۳ ۹ 
بهار ۱۷۲۰۸۱ | ۲۲,۱۲۳ ۵۲۸ ۱2/۵6 ۲۳۳۹۲۶ و۰1 0 -۲۱۳۵ 1 
تابستان ۱۳,۱۰۳ ۱ -۱۳۵۱, 16 ۱۹۱۸۸۲ ۸۶ -۷۳ 16 ۱۵۸۵۰۸ ۳۸۱۱۶۲۲ 16۳۷/۰۵ 
پائیز ۲۶,۰۱۹ | ۸۶۳۲۹ ۲۲ -۰6 ۲۲ ۱ ٩2‏ -۳۸۵ ۲۹ ۱۳ ۱۸۳۵ 
سالاثه -1۱۷۸۶۱۹ ۱۳ -1۶9۷۵ ٩۲۱۵۰۸‏ ۱3 ۱۱۸۸/۵۸ -۱۷,۱۱۵۱ ۱۹ -۱۹۵۹۵, 


جداول ۵ و نموداره نشان می‌دهند که ميانگین بارش در ایستگاه تبریز در اکثر ماه‌های دوره‌های آینده کاهش 


خواهند داشت. مقدار بارش در ماه‌های فوریه و دسامبر و در تمام طول دوره و در تمام سناریوهای اجرا شده 


سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی... ۱۷ 


سبت به دوره پایه افزایش خواهد داشت. نتایج تستان می‌دهند که میانگین بارش در ماه‌های می تا سپتامبر یعنی 
در فصول بهارتا اوسط تابستان در دوره‌های آینده با کاهش محسوس همراه است. بخصوص در دوره ۲۰7۹- 


ان کاهن توش گر اسگ 


350 350 
۹2۳45 ۹۰826 
+ 300 300 
250 250 
20 200 
مور 2 مور 3 
ی 1 
ح‌ 2 
ِ 3 
ك‌ ك‌ 
۱۱ ۲ ۱ 
هد بل رز لآ اش , بل دنه نا لاش , 
ی عم ی ی کي + کی ی دی ی ۵ ۵ ی ی ی ی کي کم ی ی ی 
ی کی شم طيي کي شیم شیم #۳ کي له کي شب يکي دی ی کی دشن کم کم شب کی کب یی اند کب اهب يم کي کی کر 
کي ی ئِ کی کم کم کل کی یی ۳ ۳ کي ده کي کي کر ۵ ی 
نم گ يب لو ۳ گبي تن کو ۲ ۵ کي کم گ ی لو کي رگن که ۵ ت ی کی 
دوره پایه ال ۲۰۹۹-۲۰۷۰ ۲۰۶۹-۲۰۴۰ 8 ۲۰۳۹-۲۰۱۰ تا دوره پایه * ۲۰۹۹-۲۰۷۰ ۲۰۶۹-۲۰۴۰ 8 ۲۰۳۹-۲۰۱۰ نا 
350 
۱۱۳۹5 
300 
250 
200 
۳ 
150 3 
حم 
« 7 
نمی 
۱ ۲ 
مد با تال له یا بط , 
۵ ین هم کی مب جي دی تم دم تم ی ی 
تم کي کم کي گم ۳ ی ات کي ۳ ند 
۳۴ سِ کي لو ۳ یر س کی تلم س ۳ ی 
0 ف 
دوره پایه 8 ۲۰۹۹-۲۰۷۰ ۷ ۲۰۶۹-۲۰۴۰ ۲۰۳۹-۲۰۱۰ 


شکل ۵ تغییرات بارش تبریز در دوره‌های ۲۰۷۰-۲۰۹۹ و ۲۰۰-۲۰۹ ۲۰۱۰-۲۰۳۹ در مقایسه با دوره پایه 
۱۹۰-۰ میلادی 

- بررسی تغییرات دمای حدافل 

در جدول + مقادیر تغییرات میانگین ماهانه حداقل دمای ایستگاه تبریز برای دوره پایه و سه دوره آینده نشان 
داده شده است. نمودار داده‌های مربوطه جهت مقایسه در شکل * آورده شده است. مطابق نمودار" و جدول ٩‏ 
تغییرات میانگین حداقل دمای ایستگاه تبریز در کلیه ماه‌ها به‌غیراز ماه نوامبر و دسامبر در دوره‌های آینده افزایش 
خواهد داشت. در دوره ۲۰۱۰-۲۰۳۹ این افزایش دما محسوس نیست ولی در دوره ۲۰۶۰-۲۰3۹ و ۲۰۹۹- 
۲۰ این افزایش کاملاٌ معنی دار و مشخص است. باتوجه به جدول مشاهده می‌کنیم که به‌طورکلی حداقل دما 


در سه سناریوی مورد بررسی برای سه دوره مورد مطالعه افزايش می‌یابد که کمترین افزایش دمای حداقل در 


۱۳/5 جغرافیا و مخاطرات محبطی شماره پیست و یکم 


دوره اول تحت سناریو ۲002.6 و بیشترین افزایش دمای حداقل در دوره وق تلاوت 35 رخ 
خواهد داد. همچنین حداقل دما به‌ طور کلی در تمام فصول افزایش می‌پابد که در فصل تابستان تا ۸ درجه یز 
افزایش دما مشاهده و کرا3ه 

حدول ۱ تغییرات دمای حداقل تبریز در دوره‌های ۲۰۷۰۹ و ۰۲۰۶۰-۲۰۹۱۹ ۲۰۱۰-۰۹ در مقایسه با دوره 


پایه ۱۹۱۰-۱۹۹۰ میلادی 


دوره 10۳2۹6 100045 1۳009۹5 
ِ کت کت ِ کت کت ِِ کت کت 
1 و ‌‌ 1 1 ‌‌ 1 و و 
سا - - - - - - - - -« 
رب 1 لِ ٍ 1 ۲ 1 1 1 1 
سن دی ح نم دی ح نم دی ح 
‌ 


ژائویه ها ۱۲۱۸۸ | ۲۶ | ۸۳۵ | ۳۵ ۵۲۱۳۵ از ۱۵ از ۵ ۸۲۱۸۲۶ 
فوریه ۳ ۳۱۹۱۸۵ از ۳۸ | ۲۸۱۳۶ ۱۱۸ از هه از ۳۸۰۵۲۵ از ۵۱۹ ۳۱ 
مارش ‏ | ۱/۱۰۸۷ از ۱۵۸۳۲۵ ۱۳۵۵4 ۲۱۳۵ | ۲۱۲۵۳۲ از ۲۸۵۸ از ۱۹۵۵۵۸ ۱ ۳,۵۰۲۵۱۷ ]۵,۹۵0 
آوریل ۲ ۸ ۲۱۰۱۳۶۵ ]1 ۱۱۷/۶۸۸۵ ۱ ۱,۵۷۷۶۲۹ ۲ ۵ ۷ 1۰۵۲ 
می هه ۱۱۲ | ۲۵ از تا ۵۲ ۰۹۸ ۰9۳۱۳۵ ۰,۹۷ ۲۹۳۹۸ 
ژون ۱۳۹/۸ ۱۳۸ 32 2-۳ ۳ ۲"( و۱۵ ۱ ۸,۲۶۵ 
حولای | ۳۵۱۹۷۹۹ | ۶۲۳۹۵۲۶ ۸ ۸ ۲۲ ۵۱۳۳۸ ۳۸۲۱۵ ۵۵۹۵۶ | ۸۳۵۳۹۳۷ 
اگوست ۱۳۶ ۱ ۱۹ ۱۳۸۹ ۳ ۱ ۸ ۱۳۹ ۷ ۱۷۱۹ 
سیتامپر | ۳,۱۰۲۱۱ ۱۲۸۹۰ ۲ ۶۶۵۱۲ ۳۳۱۳ ۸ از ۲۵ از ۵۸۲۵ ۱۸۱۱۱ 
اکتبر ۰,9۱۷ ۰۱۵۷ 2۳۳۳۹ ۸ ۱۱۳۵۸۲ ۸ ۲۱۱۳ ۵۲۸۲۵۹ 
نوامبر ۲۱۵۱۲ | ۱۱۱۲۶ | ۲۵۲۳۵ | ۱۸۰۲۸ | -۱۸۹۲۷۲ ]| ۱۱/۳۸۵ | ۲۸۵۸۷۵۵ | ۰۸۵۲ ] ۱۳۸۸۸ 
دسامیر | ۱8۱۵۷۸ ]| ۱۵۲۶۸۲ | ۲,۱۱۹ | ۱۱۱۱۷۸۸ ] ۱۱۱۲ | ۰,۵۵۵ | ۱۸۱۸۵۲ | ۰,۵۷۸ ] ۰۷۹۱۳۳۹ 
زمستان ۱۳ 4 2 ۷( ۱-۹ تک م۱۳ ۸ ا ‏ 5 
بهار ۷ ا ۱۱/۵ أ ۱۲۶ ۱۲۳۱ ۱ ۲۱۱۰۱۵۸ ۱۰۹۳۶۱ ۶۱ ۱ ۶,۱۳۸ 
تابستان | ۳۱۳۵۵۸ | ۳۸۵۸ | معا ۱۳۱۵ ال ۲۵ | ۵۵۲ | ۳۱۲ | ۵۵۶۲ ۷۸۸۲۵۵۲ 
پائیز ۰۳۳/۱ ا ۰۳۳۵۹۵ | ۰۷۸۷۵۵۳ ا ۱۵۵۳۵ ۱۳۰۸۱۲ از ۲۱۱۵۵۵۵ ۶,۱۷۵ 


سالاه |  ۱۵۵۲۵  ۱۵۲۸۸۲‏ ع ۲۱۳۳۱‏ ۲۱۵۸ ا فتهفه۲ از ۲ از ۱ ۳۵۱۱۲ ۵۲۱۸ 


سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی... ۳ 


2۹28۹5 ۱28026 


دما (درجه سانتی گراد) 
دما (درجه ساننی گراد) 


دما (درجه سانتی گراد) 


۲۰۷۰۹ مت ۲۰۴۰-۲۰۶۹ مسسلرس ۲۰۱۰-۲۰۳۹ - سقت < دوره پایه هس 


شکل 1 تغییرات دمای حداقل تبریز در دوره‌های ۲۰۷۰-۲۰۹۹ و ۲۰2۰-۲۰۳۹ ۲۰۱۰-۲۰۳۹ در مقایسه با دوره پایه 
۱۹۰-۰ میلادی 

- بررسی تغییرات حداکثر دما 

تغییرات میانگین حداکثر ماهانه ایستگاه تبریز نیز در دوره‌های پایه و آینده بر اساس سناریوهای مورد 
بررسی در جداول ۷ آورده شده است. جهت مقایسه بهتر نمودارهای ماهانه نیز آورده شده است. طبق جدول ۷ 
و نمودار ۷ میانگین دمای حداکثر دمای تبریز مثل حداقل دما در دوره‌های آینده افزایش محسوس خواهد 
داشت. در جدول ۷ مقادیر تغییرات میانگین ماهانه حداکثر دمای ایستگاه تبریز برای دوره پایه و سه دوره آینده 
نشان داده شده است. نمودار داده‌های مربوطه جهت مقایسه در شکل ۷ آورده شده است. نتایج نشان می‌دهند 
که تیال شان گنه دب کف قمع انگاه ند هر که فاها ای ناسا کر تراس‌هر رای اشاه 
افزایش داشته است. باتوجه به جدول مشاهده می‌کنیم که به‌طور کلی حداکثر دما در سه سناریوی مورد بررسی 
برای سه دوره مورد مطالعه افزايش می‌یابد که کمترین افزایش دمای حداکثر در دوره اول تحت سناریو 1004.5 


و بیشترین افزایش دمای حداکثر در دوره خر و تحت سناریو 5 رخ خواهد داد. همچنین حداکثر دسا 


۱۷۰ جغرافیا و مخاطرات محبطی شماره پیست و یکم 


به‌ طورکلی در تمام فصول به‌غیراز فصل پائیز افزايش می‌یابد که در فصل تابستان تا ۱ درجه نیز افزایش دما 
مشاهده می‌گردد. 


حدول ۷ تغییرات دمای حداکثر تبریز در دوره‌های ۲۰۷۰۹ و ۰۲۰۶۰-۲۰۹۱۹ ۲۰۱۰-۰۹ در مقایسه با دوره 


پایه ۱۹۱۰-۱۹۹۰ میلادی 


دوره 1024 ۱ ری م۱ 
۶ گ ۶ گ ۶ ۶ گ 
ستاریو | . به : 7 1 1 1 1 1 1 
نم دی ح نم دی ح نم دی حز 


ژانویه ۱ ۱۳۵۹۷ ۱ اا ا ۳۲۸۵۵ ۱۵۵ ۳۶۵۱۸۳۵ ۶1۵۱۵2۷ 
فوریه ۳۹ ۱( ۱۸۳ و۱۳ ۱۳۹ 2,9۲ ۱۳۹ ۸ ۱ ۱( 
مارس ۱۳ ۹ ۳ ۵:۲:۳۶۱ ۱۸۸ ۵( ۲ ۲ ۱۵ ۷,۵۳۱ 
آوریل سا ۱ ۳۸۵۸۹۵ | ۱۳۸۱۲ | ۳ ره | و۲ | ۳۲۵۰۱۳۱۶ | 1۵۳۶۵۳۲ 
می ا ‏ ۵ عم ۱۱۵۱۲ زا ۸ از ۱۶ ۲۵ از ۱۷۳۲۹۵۵ ۱ ۱۳۰۸۱2۵ 
ژون ۶ ۸۱۵۲۸۲ 1 ۱۰,۲۸۶۲۹ ۵ ا ۵۸۵ ۱ ۱۱۸۸۰۳۰ ۵ ۱۰۱۷۸۹۵ ۱ ۱۵,۱۸۷۳۲ 
جولای | ها م۵ ۱۳۵۵۲ ۲۳۳۸۲ ۷۱۱۱۶۲۵ ا ۵ / ( ۳ 1 ۱۱,۰۱۷۵۵۲ 
اگوست ۱۹۳ ۱۳۸۲ ات۱۳ ۱۳۵ ۱-۳ ۶( ۱۹( ی ۳ 5 
سپتامیر 0 ۰۵۳۵۲ | ۲۵۵۳۲۵ [ ۰۱۸۱۳ ۰۸۸/۳۸ ۱۱۷۹۱۷۲ ۳3۹۹ ۱۱۰۷۸۳ 
اکتبر ۱۷۵۷۹ | -۲,۲۹۵۰۱۲ | -6۷۱۸۲,ع | ۱۸۷۹۵۲ | ۲۸۲۹۰۶ | ۳۳۸۵۹۸ | ۱:۲۳ | ۰,۲۸۵ | و۰۲۰ 
نوامیر ۰۱۶۱۱ ۰۸۷۸۷۰ | ۱۵۳۲۲ | ۵۸ ]| ۰,۵۸ | -۰۲۲۹۷ ۰,۵۷۳ | ۰۲۱۸۷۲۲ | ۱۵2۰۹۵ 
دسامیر ۱,۹ ۸ ۱۰۲۹۶۶۸ ۱ 
زمستان ۱/۳۶ ۲ ۵ ۱۸۸/۰۵۵۵ ۵ ۳۳۵۹۵ ۸۲ ۳,۰۸۹ | ۶/۳۹۹۳ 
بهار ا ۳ ۵ ۱,۷۲۱ ۱ ها ۲۵۵ | ۱:۵۸ | ۵۲۲۹۵ ۵۵۱۷۷ 
تابستان ۲ ها ۶ از ۲۸۳۵ از ما ۲۵۵۵ ۱ ۵اه ها م۱۰۸۵ 
پائیز ۰۱۳۵۱۸۵ ۱۱۳۸۵۸۲ ] ۳,۵۱۸ ۰۱۱۸۵۸ ۱۱۱۷/۱۳ -۱,۳۵۵۵ و۰۳۷۵ -۰,۰۱۹۵ ۷ 


سالاه ] ۲۱۲۳۱۵۲ ۲۲۱۱ زا ۲ از ۱۵ ۳ ۳ ات۱ ۱۲ 


۵ ۵ ۰ ۰ ‌ 
سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی.. 2۸ 
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شکل ۷ تغییرات دمای حداکثر تبریز در دوره‌های ۲۰۷۰-۲۰۹۹ و ۲۰۶۰-۲۰۹۹ ۲۰۱۰-۲۰۳۹ در مقایسه با دوره 
پایه ۱۹۱۰-۱۹۹۰ میلادی 

۵ نتیجه گیری 

روش‌های مختلفی جهت ریزمقیاس کردن خروجی شبیه‌های گردش عمومی جو وجود دارد. که هر کدام از 
منطق خاصی پیروی می‌کند. از اين ری می‌توان انتظار داشت که استفاده از روش‌های مختلف نتایج متفاوتی را 
به دنبال داشته باشد که یکی از مهم‌ترین دلایل آن عدم قطعیت‌ها می‌باشد که بر نتایج نهایی تأثیر می‌گذارند. اين 
منابع عدم قطعیت عبارتند از: عدم قطعیت مربوط به نمایشنامه‌های تغییر اقلیمی عدم قطعیت مربوط به شبیه‌های 
1 عدم قطعیت روش‌های ریز مقیاس نمودن. و عدم قطعیت شبیه‌های شبیه‌سازی. 

در تحقیق حاضر اقدام به شبیه‌سازی پارامترهای اقلیمی دما و بارش با استفاده از مدل چندگانه 
خحطی ٩241۷‏ و مدل‌های گردش عمومی جو در شهر تبریز شد. در این تحقیق از خروجی مدل 620۳5۷2 
تحت سناریوهای 5,13۳4.5,6۳2.6 36۳8 برای دوره‌های آتی در قرن ۲۱ شبیه‌سازی و مورد بررسی قرار 
گرفت: 

نتایج نشان دادند که داده‌های دما همبستگی بهتری با داده‌های مشاهداتی (در قیاس با داده‌های بارندگی) 


۱۷۳ جغرافیا و مخاطرات محبطی شماره پیست و یکم 


احتمالی نرمال می‌باشد. از جمله عوامل کاهش همبستگی در بارندگی این است که عوامل مختلفی بر روی 
بارندگی موثر می‌باشد و از طرف دیگر بارندگی یک متغیر گسسته است. این نتایج با نتایج سجاد خان و 
همکاران (۲۰۰7): ساروار و همکاران (۲۰۱۰) و نوری و آلام (۲۰۱۶) همخوانی دارد؛ بنابراین حل مشکل 
همبستگی در توسعه مدل‌های تغییر اقلیم در آینده می‌بایست مد نظر قرار گیرد. همچنین نتایج این تحقیق با 
نتایج اکثر محققانی همچون کارآموز و همکاران (۲۰۰۹). سجاد خان و همکاران (۲۰۰7) ساروار و همکاران 
(۲۰۱۸) میت و همکاران (۲۰۱۲)» جیما و همکاران (۲۰۱۳) و دهقانی پیوز و هعکاران (۱۳۹۰) در آیتکه کنه 
مدل ]51251۷ توانایی مناسبی در کوچک مقیاس نمودن داده‌های دما و بارش دارد. همسو می‌باشد. تغییر اقلیم 
می‌تواند باعث تغییرات زمانی و مکانی متغیرهای اقلیمی شود. خصوصیات این متغیرها می‌تواند اثرات زیان 
باری بر روی اجزای اکوسیستم داشته باشد. با توجه به نتایم بدست آمده مشخص گردید که در طول قرن ۲۱ 
دما رو به افزایش و بارش رو به کاهش است. در ایستگاه تبریز به‌طورکلی بارش در سه سناریوی مورد بررسی 
برای دو دوره ۲۰۱۰-۲۰۳۹ و ۲۰۷۰-۲۰۹۹ کاهش و برای دوره ۲۰8۰-۲۰3۹ افزایش می‌یابد. همچنین بارش 
به‌طورکلی در فصل زمستان افزايش و بقیه فصول با کاهش بارش مواجه خواهد بود. این نتایح با نتایج گل 
محمدی و مساح بوانی (۱۳۹۰) که دوره ۲۰۶۰-۲۰۹۹ را دوره با افزایش بارش در حوضه قره سو معرفی نموده 
هم خوانی داشته ولی با نتایج رجبی و شعبانلو (۲۰۱۲) که با استفاده از مدل 51051۷1 در منطقه کرمانشاه در 
غرب ایران دورة ۲۰۷۰-۲۰۹۹ خشک‌تر دانسته‌اند. اختلاف داشتهء بطوری که بر اساس نتایج تحقیق حاضر در 
ایستگاه تبربز دوره ۲۰۱۰-۲۰۳۹ خشک‌تر از دوره‌های دیگر خواهد بود و این تغییرات بارش در 
سناریوی18)/۳4.5 بیشتر از سناریوهای دیگر می‌باشد. تغییرات میانگین حداقل دمای ایستگاه تبریز در کلیه 
ماه‌ها به‌غیراز ماه نوامبر و دسامبر در دوره‌های آینده افزایش داشته است. حداقل دما در سه سناریوی مورد 
بررسی برای سه دوره مورد بررسی افزایش می‌یابد. همچنین حداقل دما به‌طورکلی در تمام فصول افزایش 
می‌پابد که در فصل تابستان تا ۸ درجه نیز افزایش دما مشاهده می‌گردد. میانگین حداکثر دما در سه سناریوی 
مورد بررسی برای سه دوره مورد مطالعه افزايش می‌یابد که نتایج با نتایج محققانی همچون محمدی و همکاران 
(۱۳۸۹) که بر روی ایران پژوهش انجام دادند و به افزایش دما در طول سه دوره اقلیمی مورد بررسی براساس 
سناریوهای مختلف اذعان نمودنده مطابقت دارد. حداکثر دما به‌طورکلی در تمام فصول به‌غیراز فصل پائیز 
افزایش می‌یابد که در فصل تابستان تا ۱۱ درجه نیز افزایش دما مشاهده می‌گردد. همچنین هر چه به اواخر قرن 
۱ نزدیک می‌شویم این شرایط تشدید می‌شود که این امر نشان می‌دهد که وضعیت تغییر اقلیم در منطقه مورد 
نظر جدی است هر چند برای اطمینان بیشتر از بروز تغییر اقلیم در منطقه مطالعه. نیاز به مطالعات بیشتر 


تن 


سال ششم شبیه سازی تغییرات دما و بارش ایستگاه سینوپتیک تبریز طی... ۷ 


روند نزولی میانگین بارش همراه با افزیش دما باعث افزایش تبخیر و ایجاد تغییراتی مهم در کیفیت و کمیت 
منابع آب و به تبع آن تغییرات در کم و کیف محصولات کشاورزی خواهد کرد این مهم لزوم تغیبرات برنامه 
ریزی را در بهره برداری از منابع آب و کشاورزی مشخص می‌کند. برنامه‌های آینده باید به گونه‌ای باشند که 
تغییرات پیش رو اثرات زیانبار کمتری را متوجه بخش آب و کشاورزی در این منطقه از کشورمان کند. لازم 
است هم اکنون تدابیری در جهت اصلاح نظام آبیاری. روش‌های کاهش تبخیر و استفاده مجدد از فاضلاب و 
اصلاح الگوی کشت اتخاذ گردد. 
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